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Cosmological constraints from the cluster contribution to 
the power spectrum of the soft X-ray background. New 
evidence for a low erg 
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ABSTRACT 

We use the X-ray power spectrum of the ROSAT aU-sky survey in the R6 band (w 
0.9-1.3 keV) to set an upper Hmit on the galaxy cluster power spectrum. The cluster 
power spectrum is modelled with a minimum number of robust assumptions regarding 
the structure of the clusters. The power spectrum of ROSAT sets an upper limit on 
the — plane which excludes all the models with (j% above erg = O.Sri^" "^^ in 
a flat ACDM universe. We discuss the possible sources of systematic errors in our 
conclusions, mainly dominated by the assumed — T relation. Alternatively, this 
relation could be constrained by using the X-ray power spectrum, if the cosmological 
model is known. Our conclusions suggest that only models with a low value of erg 
(ug < 0.8 for rim = 0.3) may be compatible with our upper limit. We also find that 
models predicting lower luminosities in galaxy clusters are favoured. Reconciling our 
cosmological constraints with these arising by other methods might require either a 
high entropy floor or wide-spread presence of cooling flows in the low-redshift clusters. 

Key words: cosmological parameters, galaxies:clusters:general 



X 
S3 



1 INTRODUCTION 

Galaxy clusters are a useful probe of cosmological models. In 
recent years there has been a wide variety of work based on 
the study of the population of galaxy clusters as a way to 
constrain cosmological models. The mass function derived 
from N-body simulations as well as analytical approaches 
like the Press-Schechter mass function show a strong depen- 
dence of the mass function on the cosmological model. By 
using the local abundance of clusters, it is possible to set 
strong constraints on the erg — Q.m plane. 

However, both parameters are degenerate and one needs 
to go to higher redshifts in order to break this degeneracy. 
Combining the low redshift mass function with the high red- 
shift mass function it is possible to break the degeneracy in 
(78 — f^m and determine both parameters (Bahcall & Fan 
1998). Nevertheless, a measurement of the mass function 
(especially at high redshift) is very difficult to achieve since 
the mass of a cluster cannot be measured directly but re- 
quires modelling of the observations. Direct measurement of 
the mass requires the use of sophisticated techniques such 
as gravitational lensing (and high quality data) or tracing 
the mass via the X-ray emission under the assumption of 
hydrostatic equilibrium. 



An alternative approach can be devised by transform- 
ing masses into temperatures or luminosities and using the 
equivalent temperature function or luminosity function in- 
stead of the mass function. The drawback of this last ap- 
proach is that in this case a.T — M relation or a. L^ — M rela- 
tion needs to be assumed. These scaling relations suffer from 
scatter and they are not yet very well established. Even in 
the case where these relations are well known, there are still 
several systematic effects which are difficult to be accounted 
for. For instance, the luminosity of a cluster depends on the 
cosmological model since luminosities are obtained from the 
measured flux times the square of the luminosity distance. 
In the case of X-rays, the X-ray emission is concentrated in 
the centre of the cluster (due to the dependence with n 
being the electron density) and most of the emission coming 
from the external parts of the cluster is not detected. In or- 
der to compute the total luminosity of the cluster one needs 
to assume a profile for the electron density and extrapolate 
the observed emission to the outer parts. When using the 
mass, temperature, flux or luminosity functions in cosmo- 
logical studies, a precise knowledge of the selection function 
of the survey is required since the observed function will be 
compared with the theoretical one in which the integration 
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Figure 1. ROSAT all-sky survey (RASS) map in the R6 band in galactic coordinates. The map has been rcpixelized using HEALPIX. 
The colour scale is logarithmic and the units are 10~^cts/s arcmin^ . The two horizontal lines mark the two constant Galactic 
latitudes at |6| = 40°. The thick lines show the galactic meridians I = 70° (left) and I = 250° (right). 



of the mass (or temperature, flux, luminosity) starts at a 
minimum value (as a function of z) defined by the selec- 
tion function of the survey. Usually this selection function is 
inhomogeneous due, e.g., to different exposures in different 
sky directions. This makes the modelling of the theoretical 
mass (temperature, flux or luminosity) function very com- 
plicated. Also due to the selection function of the survey, 
the number of clusters observed (i.e. detected above some 
detection threshold) represents only a small fraction of the 
population of clusters. This means that cosmological studies 
based on the mass (temperature, flux of luminosity) func- 
tion are based only on the tip of the iceberg and therefore 
subject to large variations due to small number statistics. 
In this work we propose the use of the power spectrum as 
an alternative cosmological test which is less affected by the 
problems pointed out above. Our attempt is not to mea- 
sure the power spectrum of galaxy clusters (see for instance 
Schuecker et al. 2001) but to use the power spectrum of 
the diffuse X-ray background as an upper limit (see Fabian 
&: Barcons 1992 for a discussion of the X-ray background). 
This approach has several advantages over previous work. 

(i) First, since we use the diffuse power spectrum as an up- 
per limit, we do not need to detect the clusters. The only 
thing we need to do is to model the X-ray cluster power 
spectrum and compare it with the measured one. The clus- 
ter power spectrum must remain below the one due to the 
global fluctuations of the diffuse cosmic X-ray background 
(Barcons & Fabian 1988, Carrera et al. 1997, Barcons et al. 
1998, Yamamoto & Sugiyama 1998). 

(ii) In modelling the cluster power spectrum we do not need 
to know the selection function since all the clusters (detected 
and undetected) will contribute to the power spectrum. 

(iii) The previous point has a very interesting feature. We 
are not only sensitive to the tip of the iceberg but to the 
whole iceberg itself! This fact has interesting applications in 
studies like pre-lieating (see e.g Voit et al. 2002) which are 
particularly important in the low-mass regime where clus- 
ters are difficult to detect. 

(iv) Since we do not need to detect the individual clusters 
we are not affected by any extrapolation of the measured 
flux (needed to get the total luminosity or mass) making 
our results more independent of the assumed profile of the 
clusters. 

(v) The power spectrum is scale dependent. Previous work 
has used the energy spectrum as an upper limit to constrain 
the emission from galaxy clusters (see e.g Wu et al. 2001). 
However, in these studies the measured emission is contami- 
nated by the diffuse emission coming from the Galaxy. This 
is equivalent to studying just the monopole of the multipole 
decomposition. On the contrary, the power spectrum can 
disentangle the diffuse from the point-source emission since 
the diffuse emission will contribute significantly only to the 
large scales (low multipoles) while the emission from point- 
sources will dominate the smaller scales (intermediate-high 



multipoles) . 

In the next sections we will describe briefly the power 
spectrum of RASS diffuse emission (section 2), we will define 
the model to compute the cluster power spectrum (section 
3) and we will constrain the cosmological model (section 4) . 
In section 5 we revise the assumptions made in the previous 
sections and check how the results change when we change 
the model. Finally, we discuss our results in section 7 and 
summarise in section 8. Throughout the paper (unless oth- 
erwise stated), the Hubble constant is given as 100 h km s~^ 
Mpc. 



2 THE POWER SPECTRUM OF ROSAT 

The use of the spherical harmonics decomposition of the 
intensity of the X-ray background as a means to measure 
or constrain cosmological parameters was first proposed by 
Lahav, Piran & Treyer (1997). The emphasis was mostly on 
the contribution of the clustering of point sources (mostly 
AGNs). They focused their analysis on the low multipoles 
{I < 10). That sort of study has the disadvantage that the 
physics that leads to the formation and clustering of these 
exceptional sources is not simple and it is hidden in an 
unknown bias parameter. Indeed, Treyer et al. (1998) con- 
ducted that analysis on the 2-10 keV data from the HEAO- 
1 A2 experiment and found some interesting constraints on 
the bias parameter of X-ray sources, but no direct implica- 
tions on the cosmological parameters themselves. The use of 
a softer band (but not heavily contaminated by the Galaxy) 
to constrain the cluster power spectrum and derive cosmo- 
logical parameters is a much more promising task, as the 
physics of the cluster formation is much simpler and closer 
to the Cosmological model. 

The harmonic power spectrum of the diffuse ROSAT 
All-Sky Survey (RASS, see Snowden et al. 1997) X-ray back- 
ground has been computed up to Z « 1000 in a previous work 
(Sliwa et al. 2001). In that work the authors computed the 
power spectrum in different areas (subsamples) of the sky 
centred in the north pole. They also studied the effects of 
the data binning in pixels and calculated the error estimates 
due to limited counting statistics, instrumental background 
subtraction, and cosmic variance. 

On the other hand, the power spectrum of clusters in 
X-rays has been recently measured using ROSAT data (see 
Schuecker et al. 2001 for an estimation based on the RE- 
FLEX survey and Borgani et al. 1999 for an estimation of its 
Fourier transform, the two point correlation function, from 
the XBACs cluster sample). That measured cluster power 
spectrum has been also used to constrain the cosmological 
model (Schuecker et al. 2002). In that work, the modelling 
of the observed cluster power spectrum (above certain lim- 
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Figure 2. ROSAT power spectrum in tlie R6 band computed 
witli HEALPIX in different areas of the sky. Top dashed hne 
Ci, for the all-sky RASS data. Dotted hnes, Ci for the frac- 
tion of the sky with b > 40° (top) and 6 < -40° (bot- 
tom). Dot-dashed line, Ci for the part of the sky defined 
by |6| > 40° and thick solid line, C; for the optimal area: 
I e [70°, 250°] and b > 40° (equal to sample D in Sliwa et 
al. 2001). The power in this area does not show a significant 
galactic contamination. Our estimations are comparable to 
the ones obtained in Sliwa et al. (2001). 



iting flux) is a complicated procedure which involves several 
assumptions (such as the — M relation which was taken 

from Reiprich & Bohringor 2002). Results obtained in a dif- 
ferent study (Viana et al. 2002) seem to indicate that there 
could be an inconsistency between the REFLEX luminosity 
function and the Lx — M relation obtained by Reiprich & 
Bohringer (2002). 

In this work, instead of detecting clusters and finding the 
power spectrum from the sample of detected clusters, we 
will use the RASS power spectrum to set an upper limit 
on the cosmological model. This approach has several ad- 
vantages outlined in the previous section. In brief, we recall 
that the main advantage of this work will be a better control 
of the systematic errors. 

We have re-estimated the power spectrum in an indepen- 
dent waj' (compared with the estimation of Sliwa et al. 
2001). In our process, we have rcpixelized the RASS data 
using the package HEALPIX* and computed the power spec- 
trum with the same package. HEALPIX was originally de- 
signed to analyse mega-pixel CMB data and is optimised 
for the computation of the harmonic power spectrum. 
The pixel size of the RASS maps is 12 arcmin and the maps 
are given in different bands. We choose the band R6 which 
is the one in which the cluster signal is expected to be at its 
maximum compared to the other components (local bub- 
ble, galactic diffuse emission, extragalactic AGN's). This 
band is also the best in terms of the instrumental response 
and background contamination. We repixelise the RASS in 
HEALPIX using a pixel size of « 13.7 arcmin (Nside = 256) 
and compute the power spectrum with the HEALPIX sub- 

* available at http://www.eso.org/science/healpix. Copyright 
1997 by Eric Hivon and Krzysztof M. Gorski. All rights reserved. 



routine anafast. Nside = 256 gives the closest HEALPIX pix- 
elisation to the original ROSAT's one. Due to the large pixel 
size, we can compute the power spectrum up to only I = 500. 
Beyond that point we do not have information about smaller 
scales. With this pixel size, the effect of the PSF of ROSAT 
becomes negligible and we will not include it in our analysis. 
However, at smaller scales, the ROSAT PSF can have an 
important effect on the measured power spectrum. 
In figure 1 we show the RASS data in the R6 band af- 
ter repixelisation with HEALPIX. Several features can be 
observed in the map, the most prominent one being the 
Galactic plane. The diffuse X-ray maps used in this work 
have been cleaned of the most prominent point sources (see 
Snowden et al. 1997). The brightest point sources (Voges 
et al. 1999, Voges et al. 2000) were removed making use of 
the full resolution of the PSPC detector (1.6' x 1.6' pixel 
size) using a minimum source-excision radius of 3.5'. 50000 
sources where removed above a threshold of 0.02 cts/s in 
the R6-I-R7 band. Below this threshold, there are still many 
point sources which could contribute to the power spectrum. 
In the maps, there are are several strips with no data (black 
bands in figure 1). There are also many pixels with count- 
rate, specially in the southern hemisphere. 

We arc interested in computing an estimate of the power 
spectrum in an area of the sky which is large enough to avoid 
cosmic variance fluctuations but also that has a low fore- 
ground contamination and low noise level. The contamina- 
tion level is determined basically by the amount of Galactic 
emission. We should then avoid regions with strong galac- 
tic emission. The noise level is determined by the exposure 
time. We should then concentrate on those regions having 
low galactic emission and large integration times. We have 
selected a region with b > 40° and I £ [70°, 250°] (totalling 
Ri 2400 deg^) as the optimal region for our analysis (see fig- 
ure 1). Several groups have selected also this region for their 
studies of the X-ray background (Soltan et al. 1996, Miyaji 
et al. 1996, Kneissl et al. 1996, Sliwa et al. 2001). 

In figure 2 we show the power spectrum in different ar- 
eas of the sky. We have represented the power spectrum in 
different situations where we exclude (or include) the galac- 
tic plane (|6| above or below 40°). The galaxy contributes 
on large scales (low I's) but also introduces power at inter- 
mediate I's. The power spectrum calculated in our optimal 
area is shown as a thick solid line. In this case, no structure 
due to the Galaxy is seen, suggesting that this area may be 
dominated by the cosmological signal (galaxy clusters and 
faint AGNs) . The optimal area has also a very low Nh col- 
umn density, typically ~ 2 x 10^" cm~^. The absorption of 
the X-ray fiux by this column density is only a few percent 
in most of the pixels in this area, an effect that can be safely 
neglected. 



3 THE CLUSTER POWER SPECTRUM 

We have used the halo model to compute the power spec- 
trum due to galaxy clusters (Cole & Kaiser 1988, Bartlett & 
Silk 1994, Komatsu & Kitayama 1999, White 2001, Cooray 
2002, Komatsu & Seljak 2002). In this model, the power 
spectrum is due to two main contributions, the two- halo con- 
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tribution and the single-halo contribution. The first one in- 
cludes the correlation between clusters and only contributes 
significantly in the very large scale regime. The single-halo 
contribution accounts for the individual contribution of each 
cluster to the power spectrum and dominates at smaller 
scales. 

The large scales (low /) of the power spectrum in the opti- 
mal area may be affected by window effects and cosmic vari- 
ance (our optimal area is just « 9% of the sky). Although 

the power spectrum in the optimal area seems to be free 
of galactic contamination, it may happen that the galaxy is 
still contributing in this area. For all these reasons, the low 
I regime should be excluded in our cosmological analysis. 
We will concentrate only on the intermediate-small scales 
for which one can consider only the single-halo contribution 
and neglect the two-halo contribution. 
However, we should keep in mind that, as noted in Komatsu 
& Kitayama (1999), the clustering contribution can amount 
to 20% — 30% of the single halo-contribution at degree an- 
gular scales (I ~ 100) (see also Lahav et al. 1997 for smaller 
multipoles). Therefore, by neglecting the two-halo contribu- 
tion we are being conservative in our assumptions. The real 
power spectrum of galaxy clusters will be about 20% — 30% 
larger than the one considered in this work. 

The single-halo contribution is just an integrated effect 
over redshift of the individual contribution from each cluster. 



C, 



dz 



dV{z) 
dz 



,^,dN(M,z) 



(1) 



where dV{z)ldz is the volume element, dN{M, z)/dM is the 
cluster mass function and pi[M,z) is the power spectrum 
(multipole decomposition) of the X-ray 2D profile of the 
cluster with mass M at redshift z. In this work we have 
assumed the Press-Schechter approach for the mass function 
(Press & Schechter 1974) and we use a numerical fit to the 
multipole decomposition of the X-ray 2D profile as we will 
describe in the next subsection. 



3.1 The multipole decomposition of the X-ray 2D 
profile 

Since the power spectrum of galaxy clusters is proportional 
to the square of the density profile, it is important to be 
careful in the selection of this profile and make the safest 
assumptions about it. The usual approach is to consider 
a numerical model for the electron density profile (a (3- 
model or a numerical fit to N-body simulations) and com- 
pute the flux from the integrated signal of the profile. 
Then, assume isothermality of the plasma and some rela- 
tion between the mass of the cluster and its temperature. 
In this process there are several dangerous steps. Assum- 
ing a bremsstrahlung spectrum, the monochromatic volume 
emissivity of the intra-clustcr medium scales as oc r? l\fT , 
where n is the electron density and T its temperature. Then 
the integrated bolometric luminosity scales as the square of 
the central density times the square root of the temperature 
(assuming the cluster is isothermal). The central density has 
been poorly estimated so far in galaxy clusters and there is 
a lot of uncertainty in its typical value. If the cluster is as- 
sumed to be isothermal and the scaling relation between 
mass and temperature is taken to be the one predicted by 
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Figure 3. Multipole decomposition of tlie 2D X-ray profile for 
different core radii (rc). I^om left to right, = 27,13.7,7,3.5 
arcmin. Tlie dashed line represents the real multipole decompo- 
sition and the dotted line is the numerical fit of equation 5. 



the spherical collapse model, one ends up with a, — T 

relation which is inconsistent with recent observations. We 
can conclude that, modelling of the X-ray emission in galaxy 
clusters from numerical models can be a quite unsafe pro- 
cess. It is more useful to formulate the model in terms of ob- 
servational quantities rather than theoretical models which 
fail in explaining the observations. 

Basically what we need for our model is just an expression 
of the form; 



pi{M,z)=p„iM,z)*f{l,M,z) 



(2) 



where Po is the normalisation and f{l, M, z) contains the 
angular dependence of the multipole decomposition. To un- 
derstand why this simple form will suffice, it is useful to 

think in terms of a single point source. This source will have 
a constant power spectrum {f{l, M, z) = 1) and a normali- 
sation which is equal to; 



Po{M,z) = 4-n-\Mean\ 



■ 4-K 



St ' 
4-K 



(3) 



where Mean is the total signal, St, divided by the area 
of the sky. In the case of a cluster, the situation is similar 
except that in this case the source will be extended and 
f{l, M, z) will not be constant anymore. In this case, we have 
to assume a profile for the electron gas. We assumed a 13- 
model (with /3 = 2/3) profile for the electron cluster density. 
This assumption is supported by a wealth of observational 
data. We also have truncated the profile at the virial radius 
defined as = pre with p a fixed parameter (p w 10) and 
rc the core radius. 
rio 



n{r) = 



1 + [r/rc 



(4) 



Instead of fixing the central density, rio, we will fix the total 
luminosity so the value of Uo will be irrelevant in our model. 
The only free parameter is the core radius, r^- This will be 
an advantage compared to other works since we can make 
our model consistent with real observations rather than with 
simulations. 

To compile an expression for the multipole dependence of the 
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2D profile, we have computed f{l, M, z) for different clusters 
changing the only free parameter, Vc, and then fitting the 
result to an analytical form. The best fit we found is; 

f{l,M,z) = i (exp(-Ci,.J +exp(-^/a^)) 

with, 

t — ;2„1.5/(0.815+0.35r<,) 

S!,T-c ' ' C 

where the core radius, rc, is given in rads (see figure 3). 
The mass and redshift dependence of /(/, M, z) is in Tc- 



(5) 
(6) 



P 



(7) 



Although we have Eissumed the expected redshift depen- 
dence of the self-similar spherical collapse model with firn, = 
1, we will adopt this form for simplicity. As we will show 
below, the power spectrum of galaxy clusters will be domi- 
nated by the low redshift population for which the specific 
form of the redshift dependence is not very important. The 
value of ro must be assumed and we should keep in mind 
that typical core radii for massive clusters arc around 100 
kpc. We will adopt a value of ro/p = 130 h^^ kpc for a mass 
of M = \Q^^h~^MQ at redshift and will discuss other val- 
ues later. 

Regarding the normalisation of the profile, po, we only have 

to determine Mean, that is, the total signal divided by the 
area of the sky. To avoid systematic effects due to unrealistic 
modclliug based on simulations, we have parametrised the 
total signal in terms of observed quantities, like the — T 
relation. We used the scaling relations derived in Diego et 
al. (2001). In that work, the authors made a consistent fit to 
the mass function, temperature function. X-ray fiux function 
and X-ray luminosity function. In the fit both, the cosmolog- 
ical model, and the cluster scaling relations were considered 
as free parameters and were fitted to the data. Using the 
Lx — T relation we can obtain the total signal (flux) of the 
cluster, St, from its temperature. 



fc 



(8) 



'4nDi{z)^ ' 4nDi{z)^ 

where we use the values of Lo, a and tp obtained in Diego 
et aL (2001) {Lo = 1.12 x 10*^h-^erg/s, a = 3.2 and 
■i/) = 1 for a flat ACDM universe). The temperature is ob- 
tained using the T — M relation found in Diego et al. (2001) 
{Tkev = 9.48M}'5^^(1 -I- z)). This relation is different from 
the best fitting function to the observed T — M relation (see 
e.g in Nevalainen et al. 2000). However, we will use this re- 
lation to be consistent with the above — T relation and 
the constraints in Diego et al. (2001). In section 6 we will 
discuss other (possibly more realistic) alternatives. 
We should note that by changing the parameters Lo, a and -ip 
we can include also the uncertainty in the T—M relation into 
an uncertainty in the L^—T relation. That is, the parameters 
Lo, a and will somehow mimic variations in the param- 
eters of the T — M relation. This choice is not arbitrary. 
We could easily model the total luminosity as a function of 
mass and use the Lx — M relation as our free-parameter scal- 
ing relation. However, wc prefer to parametrise our model 
in terms of the Lx — T relation since it is observationally 
better constrained than the Lx — M relation. 
The only free parameters in our model are then the L^ — T 
relation and the rc modelling. We will study their effect 




Figure 4. Dependence of the X-ray power spectrum on the cos- 
mological parameters erg, Q, and V. The solid Une is the reference 
model ((78 = 0.8, = 0.3 and V = 0.2). We show the effect 
of changing one of the three parameters while keeping unchanged 
the other two. The effect of changing as while O and F are fixed to 
the values of the reference model is shown in the dotted lines. Bot- 
tom dotted line is for as = 0.7 and top dotted line is for as = 0.9. 
Dashed lines show the effect of changing in 0.1 units. Bot- 
tom dashed line is for = 0.2 while the upper dashed line is 
for f2m = 0.4. Finally, the dot-dashed lines show the effect of 
changing the shape parameter, F. Bottom dot-dashed line is for 
F = 0.25 and upper dot-dashed line for F = 0.15 



later. 

The factor /c is the conversion factor between ergcm"^ s~^ 
to 10~®cts/s arcmin'^ which are the units in which the 
RASS data are given. We compute this factor with the X- 
ray package XSPEC (Arnaud 1996). We have assumed a 
z-bremsstrahlung model and considered only the ROSAT 
PSPC channels of the R6 band (channels 91-131, see Snow- 
den et al. 1997). The ROSAT response was modelled us- 
ing the response matrix file of the detector PSPC-C. Us- 
ing a Raymond-Smith for the model spectrum with 0.3 so- 
lar metallicity in the R6 band does not introduce a signifi- 
cant difference compared to the z-bremsstrahlung. The K- 
corrcction was also included in the factor fc. This factor was 
computed for each temperature and redshift. 
As an example, XSPEC predicts a conversion of 1 cts/s = 
8.6 X IQ-^^erg/s cm-^ in the R6 band (for kT = 7 keV and 
z = 0.5) 



3.2 Cosmological dependence of C'l 

The power spectrum shows an important dependence with 
the cosmological parameters as, ^ and P. In figure 4 we 
show some examples where wc vary the normalisation of the 
power spectrum, as, the matter density, Q, and the shape 
parameter, P. The dependence with as and Q is very im- 
portant. A change of 0.2 units in leads to a change in the 
power of one order of magnitude. Something similar happens 
with as although in this case the change is slightly weaker. 
The dependence on P is not so important. In particular, in 
our region of interest {I m 100), P does not have a very sig- 
nificant effect on the power spectrum. For this reason we 



© 0000 RAS, MNRAS 000, 1-13 



6 Diego et al. 




1 10 100 1000 10000 



Figure 5. Power spectrum of galaxy clusters (with M G [3 X 
10", I0^^]h-^MQ) split in diflferent redshift intervals (erg = 0.8, 
Urn = 0.3). Solid line z e [0.01,2], dotted line z g [0.01,0.05], 
dashed line z e [0.05,0.2] and dot-dashed line z e [0.2,2]. 



will fix r in our analysis to its most favoured value T = 0.2 
(from 2dF, Percival et al. 2001) and we will vary only as 
and flm- 

3.3 Mciss and redshift dependence of C; 

Before comparing the cluster power spectrum with the 
ROSAT data in the R6 band, it is interesting to look in 
more detail at the way the different cluster populations con- 
tribute to the power spectrum. Given a cosmological model, 
the cluster mass function is just a function of two variables, 
mass and redshift. We can now examine the cluster pop- 
ulation at different mass and redshift intervals to find out 
which regimes dominate at different scales. 
Low redshift clusters will dominate the large scales since at 
low redshift the clusters appear larger (figure 5). Also, it is 
difficult to find large clusters at high redshift. They usually 
appear at low z. On the contrary, high redshift clusters will 
dominate the high-Z regime. At high z, clusters will appear 
with small angular scales. Also, the typical mass of a cluster 
decreases with redshift. However, clusters beyond z « 0.2 
will contribute only at the very small scales {I > 2000 or 
scales smaller than ~ 5arcmin). At scales larger than 5 ar- 
cmin, only the clusters below z w 0.2 are relevant in the 
power spectrum. We also show the contribution in the range 
z G [0.01, 0.05]. Many Abell clusters will contribute to this 
part of the power spectrum (min ZAbeii ~ 0.01). 

In terms of the masses (figure 6), the power spectrum 
is dominated by intermediate mass clusters (dashed line). 
The very massive clusters are rare and they do not domi- 
nate the power spectrum at any scale. They are big clusters 
and usually appear at low redshift. As a consequence, they 
will contribute significantly only to the large scales. On the 
contrary, intermediate-mass clusters are common in a wider 
range of redshifts and their integrated signal dominates the 
signal of the very massive clusters. The intermediate-mass 
clusters are bright enough to dominate the power spectrum 
at practically all scales. Finally, the low-mass clusters will 
contribute only to the very small scales (high l-s) . They are 
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Figure 6. Power spectrum of galaxy clusters (with z G [0.01,2]) 
split in different mass intervals (cts = 0.8, ilm = 0.3). Solid line 
M e [3x 10i3, lOi6]/i-iM0. Dotted line shows the contribution of 
the dark clusters, M e [3 x 10^^, 10^*]ft~^Mo. Dashed Une shows 
the dominant contribution of the intermediate-mass clusters, 
M e [10", 10^^]h-'^MQ and the dot-dashed line shows the con- 
tribution of the very massive clusters, M £ [10^^, 10^^]h~^MQ. 



more abundant than the intermediate-mass clusters but the 
signal is weaker. 

Combining the results of figures 5 and 6 we can con- 
clude that the cluster power spectrum up to Z « 500 will 
be dominated by the cluster population with z < 0.2 and 
intermediate masses, M e [lO",lO^'^]/i-^M0. This can be 
seen more clearly in figure 7 where we present the differen- 
tial contribution to different multipoles as a function of mass 
(lower set of curves) and redshift (upper set of curves). For 
representation purposes, all the power spectra has been nor- 
malised to 7 for the redshift case and to 4 for the mass case. 



4 CONSTRAINTS ON as - Qm 

After computing the power spectrum from eqn. 1 for hun- 
dreds of models where we change as and il, we can com- 
pare the predicted power spectrum with the observed RASS 
power spectrum in the R6 band and exclude all the models 
for which we have more predicted power than the observed 
one. We show the result in figure 8. The dotted line is a 
numerical fit to the upper limit. 

as = 0.50-°-'' (9) 

Our limit excludes regions in the as — plane which 
are within the confidence region of many other methods. 
In figure 9 we compare the RASS power spectrum with 
three models above, on, and below the upper limit curve 
{as = 0.75,0.8,0.85, for O = 0.3). The range of Vs which 
exclude the models above the upper limit are up to Z « 200 
which corresponds to angular scales 9 « 1°. At these scales 
only the nearby clusters (like the Abell clusters) contribute 
to the power spectrum (see figure 5 and 6). Hence, with this 
data we are not very sensitive to the population of clusters at 
moderate and high- redshift. This fact is in agreement with 
the measured cross-correlation between the angular position 
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Figure 7. Differential contribution to different multipolcs as a 
function of redsliift (top curves) and mass (thick bottom curves). 
Solid line is the contribution to the multipole I = 100, dashed 
line for I = 1000 and dot-dashed for I = 10000. All the curves 
have been normalised to 7 (top) and 4 (bottom). The numbers 
on top of the curves indicate the multipole to which they are 
contributing. The mass scale is shown in the upper x-axis and 
the x-scale in the bottom x-axis. 




Figure 8. Upper limit constraints in the as — ilm plane (solid 
line). The dotted line corresponds to the curve, as = 0.5Q^'^^ . 



of Abell clusters and the intensity of the X-ray background 
(Soltan ct al. 1996). In that work, the authors found that 
the cross-correlation extends up to several degrees (i.e. the 
cluster-X-ray background cross-correlation spreads over the 
range I < 100). 

The RASS data still contains some contribution of non- 
rcmovod point sources which could be contaminating the 
power spectrum at sub-degree scales. This fact suggests that 
a study carried out on a smaller region of the sky with a 
more careful point source removal could render better con- 



Figure 9. ROSAT power spectrum in the R6 band (solid line) 
compared with three models above and below the upper limit 
(dotted lines). Prom top to bottom, as = 0.7, as = 0.8 and 
<T8 = 0.9. All models are for Q.m = 0.3. The dashed line part of 
the ROSAT power spectrum was not used in the calculation of 
the upper limit in figure 8. 



straints on the cosmological model. This opens an exciting 

possibility for future work carried out with Chandra and/or 
especially with XMM-Ncwton data. XMM-Newton has a 
higher sensitivity for extended sources and wider field of 
view than Chandra. This makes XMM an ideal instrument 
for wide field X-rays surveys (like the XMM-LSS, Refregier 
ct al. 2002a) . However, the analysis of the arcmin-subarcmin 
scales is beyond the scope of this paper and will be discussed 
in a subsequent paper. 

The constraints derived in this section are robust since 
they were obtained under a minimum number of assump- 
tions. 

However, variations in these assumptions (Lx — T, core ra- 
dius) can have an important effect on the resulting power 
spectrum. It is important to study how the upper limits 

obtained in this section can be affected by varying our as- 
sumptions {Lx — T relation and core radius, rc). Although 
the model we assumed for — T and core radius was con- 
sistent with several optical and X-ray cluster data sets (as 
shown in Diego et al. 2001), there are some uncertainties in 
the values of the parameters assumed in the model (Lo, a, 
tp and Tc). A different set of values for these parameters will 
render a different upper limit on the as — plane. Our 
results can be therefore affected by systematic errors due to 
the wrong election of the parameters of the model. In the 
next section we will attempt to see how much does our up- 
per limit change when we make different assumptions in our 
model. 



5 SYSTEMATIC EFFECTS 

The main source of possible systematic error is the normal- 
isation of the assumed — T relation, Lo (eqn. 8). There 
is still an important uncertainty in this relation and even a 
break could be present at low temperatures. This break in 



© 0000 RAS, MNRAS 000, 1-13 



8 Diego et al. 




1 10 100 1000 10000 



Figure 10. Systematic effects. Effect of changing Lo, a and ip in 
the Lx—T relation (eqn. 8) for a (Tg = 0.8, Qm = 0.3 model. The 
thick solid line is the reference model. (Lo = 1.12 X 10'*'^/i~'^erg/s, 
a = 3.2, jp = 1). The dashed line correspond to the same value of 
Lo but a = 2.7. The dotted lines are for Lo = 2.25x10'^'^ h-'^erg/s 
with a = 3.2 (top) and Lo = 0.56 x 10*"^ h'"^ erg / s with a = 3.2. 
In the dot-dashed line we fix Lo = 1.12 x h~^10'^'^erg/s, and 
a = 3.2 and we change ip in four units, tfj = 5. 



the scaling has been interpreted as due to preheating phe- 
nomena (e.g. Volt & Bryan 2001, Babul et al. 2002). Since 
the power spectrum of the clusters is proportional to the 
square of the total signal (eqn. 3). The power spectrum will 
be therefore very sensitive to the normalisation of the L^ — T 
relation (Lo in eqn. 8) . This is illustrated in figure 10 where 
we change Lo by a factor of two and the power spectrum 
changes by a factor of four (dotted line) as expected. The 
dependence of the power spectrum on the scaling exponent, 
a, is not as dramatic as in the case of Lo. When wc change 
Q between a = 3.2 and a = 2.7 the power spectrum changes 
from the solid line to the dashed line. Finally, the factor 
controlling the redshift dependence, ip, is not very relevant 
up to Z « 100 but at smaller scales {I > 100) it can become 
important since as we discussed in section 3, at the smallest 
scales the power spectrum is dominated by the high redshift 
population. The previous plot demonstrates that one must 
be careful in choosing the L^ — T relation since the results 
can change significantly. 

The second assumption on our model is about the 2D 
profile. In this case, our only free parameter is the scaling 
of the core radius, rc. Wc will study the effect of chang- 
ing the core radius by changing the ratio Vo/p in equation 
7. The effect of changing the scaling of core radius is only 
important at small scales (large I). This is not surprising 
since at low I's the value of power spectrum is driven only 
by its normalisation {po in equation 2). When wc change the 
parameter p (ratio between the virial and core radius), we 
observe that clusters with larger p (i.e smaller core radius 
when To is fixed) have more power at smaller scales (figure 
11). This point is interesting since it shows how our con- 
straints on (78 — arc not very sensitive to the cluster's 
profile. The Lx—T relation will be the main source of un- 
certainty. 
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Figure 11. Systematic effects. Effect of changing the parameter 
p in equation 7 while fixing ro = 130h~^ Mpc. The solid line is 
the reference model with p = 10 (as = 0.77, Qm = 0.3). The two 
dashed lines are for p = 5 (bottom) p = 20 (top). As expected, 
clusters which are more extended (p = 5) have less power at 
smaller scales but do not change the power in the large scales (see 
text). On the other hand, clusters which are more clumped than 
the p = 10 models help to increase the power at smaller scales 
(without affecting the large scales) and they can produce a good 
match to the shape of the observed X-ray power spectrum (dotted 
line) at high multipolcs (0 < 30 arcmin). This fact suggests that 
the dumpiness parameter, p, could be in the range 10 < p < 20. 



From the previous results wc can conclude that due to 
our uncertainty in the Lx — T relation and the cluster density 
profile (parametrised as a function of rc) wc have an uncer- 
tainty in the theoretical power spectrum of clusters of w 1 
order of magnitude. We will incorporate this uncertainty in 
our upper limit in the us — fim plane by rescaling our theoret- 
ical model by a factor of 5 in both directions (multiplication 
and division). By doing this we will have an optimistic and a 
pessimistic model. The optimistic model will contain 5 times 
more power than the model assumed in section 4 while the 
pessimistic one will contain 5 times less power. The differ- 
ence in power between the optimistic and the pessimistic is 
then a factor 25. The result is shown in figure 12. Even in the 
pessimistic case, the constraints are still important (models 
like (Tg = 1, ilm = 0.3 can be excluded). As an example, the 
pessimistic case exclude more or less half of the confidence 
region obtained from the REFLEX cluster power spectrum 
(see figure 4 in Schuecker et al. 2002). 

It is more interesting to see what happens in the optimistic 
case. Here, models with erg — 0.8, Qm ~ 0.3 would be ex- 
cluded. This particular model is actually one of the best in 
fitting several recent data sets (CMB, 2dF, etc). This means 
that, if the underlying cosmological model is in fact erg = 0.8, 
fim = 0.3, our optimistic case is too optimistic! This opens 
the possibility of using the diffuse X-ray background as a 
way to study cluster physics rather than constraining the 
cosmological model. In the hypothetical scenario whore the 
cosmological model is known, one can use the X-ray power 
spectrum to constrain the physical parameters of the intra- 
cluster plasma. If, for the known cosmology, one observes 
that the theoretical cluster power spectrum exceeds the ob- 
served one, this means that the assumptions made on the 
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Figure 12. Upper limits constraint in the as — Clm plane for the 
optimistic (bottom solid line) and the pessimistic (top solid line) 
models. The dotted line is the numerical fit obtained in section 4 



{as = o.5n:^°-^^). 



modelling of the cluster power spectrum were wrong (too 
optimistic case). One should then reintroduce changes in 
the model in order not to exceed the observed X-ray power 
spectrum (pre-heating, electron density profile, evolution of 
the typical luminosity with redshift, baryon fraction, etc) 



6 BEYOND THE COSMOLOGICAL MODEL 

The main difference between the constraints obtained in this 
work and those obtained by other methods (based on galaxy 
clusters) is that we can have a good control on the systematic 
errors. Our results are not afi^ected by a changing selection 
function in the sky, or extrapolation of the observed central 
cluster 2D profile beyond the noise level. 

There are only two key assumptions in our model, the 
Lx — T relation and the cluster density profile. At the scales 
relevant in this work {I < 100) the uncertainty in the cluster 
density profile is not very important. On the contrary, the 
main uncertainty comes from our partial knowledge of the 
Lx—T relation. Most of the cosmological methods based 
on galaxy clusters data obtain the best fitting cosmological 
model assuming that the physics of clusters are well known. 
If the uncertainties in the cluster physics are not included, 
then a biased estimate of the best cosmological model may 
emerge as a result. We are approaching an era where the best 
cosmological model is quickly converging to a set of preferred 
values. One can then think the other way around, assuming 
we know the cosmological model, what can we infer about 
the physics in clusters ? 
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Figure 13. Lx — T relation compared with the prc-hcatcd model 
(solid line) in Babul et al. (2002). The dashed line represents 
the main model used in our analysis (figure 8). That model was 
the best fitting model found in Diego et al. (2001) (Lo = 1.12 x 
10^'^h~^erg/s, and a = 3.2, see eqn. 8) after fitting the cluster 
mass function, temperature function. X-ray luminosity function 
and X-ray flux function. Also shown are two of the models used 
in figure 10. Dotted line is for Lo = 1.12 x 10'^^h~^erg/s, and 
a = 2.7, and dot-dashed line is for Lo = 2.25 X h^^erg / s, and 
a = 3.2. These two models also produced a reasonable good fit in 
Diego et al. (2001). In he above plot we have assumed h = 0.75 
in Lo which was the value adopted in Babul et al. (2002). All 
the models represented in this plot give a reasonable fit to the 
observed Lx—T relation (see figure 5 in Babul et al. 2002). 



6.1 Preheating and the Lx—T relation 

In figure 13 we present some of the estimates of the Lx — T 
relation used in this work compared to a model including 
pro-heating (solid line, Babul et al. 2002). All the models 
represented in the plot make a good fit to the most re- 
cent estimations of the Lx — T relation (see e.g, Marke- 
vitch 1998, Allen & Fabian 1998, Helsdon & Ponman 2000 
). The pre-heating model corresponds to an entropy floor 
of kTn~^^^ ^ 427 keV cm . This pro-heated model gives an 
excellent match to the observations (see Babul et al. 2002). 
The model used in section 4 corresponds to the dashed line. 
The dotted line and dot-dashed line where used in section 5 
to illustrate the change in the cluster power spectrum when 
the parameters Lo and a where changed. The main differ- 
ence between the pre-heated model and the single scaling 
relation is that at temperatures T w 1 keV the pre-heated 
models predict a larger luminosity than the standard single 
scaling relations. Clusters with T of a few keV correspond to 
intermediate masses. We have seen that this range of masses 
dominate the power spectrum so we should expect some kind 
of dependence on the power spectrum with the amount of 
entropy in the pre-hcatcd model. As shown in Babul et al 
(2002), lower entropy floors will produce a larger shift of the 
relation to higher X-ray luminosities (lower a) while higher 
entropy floors will result in a steeper Lx — T relation (higher 
a). At kT below 1 keV the situation inverts and the pre- 
heated model predicts clusters to bo less luminous than in 
the single scaling case. Previous studies based on the inten- 
sity of the soft X-ray background suggest that there must 
be some level of pre-heating on galaxy clusters, for example 
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Figure 14. Change in power when prc-hcating is considered. The 
solid line is the reference model (eg = 0.8, Qm = 0.3, Lo = 1.12 X 
10'^'^h~^erg/s, a = 3.2 and i/i = 0.0. In the reference model the 
temperature of the cluster is computed from the T — M relation 
obtained in Diego et al. (2001) {T{keV) = 9ASM°^^{l + z)). The 
dashed line shows the power spectrum when the X-ray luminosity 
is computed from the pre- heating model (Babul et al. 2002) shown 
in figure 13 and the temperature is obtained from a numerical 
fit to the T - M relation (Ncvalaincn ct al. 2000), T{keV) » 
S.OMj'j^*". For reference (dotted curve), wc also show the cluster 
X-ray power spectrum when wc use the Ncvalaincn ct al. (2000) 
T — M relation but the non-preheated model is used to get the 
X-ray luminosity (i.e. Lo = 1.12 X W^'^h~'^erg / s and a = 3.2 for 
all T) 



from AGN heating of the intracluster gas. Otherwise, their 
contribution to the soft to X-ray background exceeds the 
derived upper hmit (see e.g. Wu et al. 2001). In figure 14 we 
show how the preheating can change the power spectrum of 
clusters. The solid line is the non-preheated model used in 
section 4 but with non-evolving — T {ip = 0) (dashed line 
in figure 13). When we use the pre- heated model of Babul et 
al. (2002) (solid line in figure 13) we see an increase in the 
power spectrum by a factor of ~ 2 (dashed line) . In the pre- 
heated model case we also have computed the temperature 
of the cluster (from its mass) using a numerical fit to the ob- 
served T- M relation (Nevalainen et al. 2000). That fit dif- 
fers from the best-fitting T—M function found in Diego et al. 
(2001). The Nevalainen's T— A4 relation is in perfect agree- 
ment with the latest estimate of the T — M scaling (Shimizu 
et al. 2002) where the authors find feT(keV) = 7.7M?i^^ 
(Mi5 is expressed in lO"'^^/i^"'^A/0). It is also in agreement 
with previous results which suggest that the T—M scal- 
ing relation differ from the expected self-similar behaviour 
(Mohr et al. 1999). When we use Nevalainen's T - M (but 
keeping the rest of parameters of the model corresponding to 
the solid line in figure 14) the new power spectrum changes 
from the solid line to the dotted line. 

If we assume that the pre-heating model used above is a 
good description of the L^ — T relation, we can extract some 
surprising conclusions. Our upper limits exclude many of the 
previous constraints on this parameter (sec Rcfrcgier et al. 
2002b for a recent compilation of constraints on erg). Prom 
the results compiled in Refregier et al. (2002b), only the 
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Figure 15. ug — ilm constraints assuming the pre-heating model 
of Babul et al. (2002) and the T — M relation of Nevalainen et al. 
(2000), (thick solid line). This upper limit would exclude models 
like IT8 = 0.8, Qm = 0.3. For reference we also show the result 
obtained in section 4 (dotted line). 



recent 2dF+CMB (Lahav et al. 2001) and galaxy clusters 
(Seljak 2002) constraints are well inside our upper limits. 
This can have two interpretations. The first interpretation 
is that those constraints derived using other methods which 
are in contradiction with the RASS upper limit on erg need 
to be revised with respect to the assumptions made in the 
appropriate analyses. The second interpretation is that we 
are being too optimistic in our assumptions. The assumed 
Lx — T relation is not a good description of the underly- 
ing Lx — T relation in galaxy clusters and we may need to 
inject more entropy in the clusters in order to make them 
less luminous (or any other mechanism which decreases the 
X-ray luminosity). This second interpretation is interesting 
because it shows us how we can study cluster physics if wc 
know something about the cosmological model. If we assume 
that the underlying cosmological model can be described by 
the parameters as = 0.8 and fi^ = 0.3 (in a flat ACDM 
universe with F = 0.2) then we can ask whether or not the 
assumed physics (Lx — T relation) is in agreement with the 
observations. We can test this by imposing that the theoret- 
ical power spectrum must be below the observed one. 

6.2 L^-T and T-M from ag - fim 

Although there is still a very lively debate about the value 
of as, there is a general consensus in the allowed range for 
as, as G (0.7— 1.1) for Qrn = 0.3 with a most favoured value 
of as w 0.8. The fact that the preheating model consid- 
ered above is excluding this most favoured value can have a 
simple interpretation, which is that our assumptions about 
the luminosities of the clusters are too optimistic. The clus- 
ters must be less luminous than we assumed in order to 
allow the upper limit to stay above the preferred value of 
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Figure 16. Upper limits on the value of Lo as a function of a for 
a cosmological model with as = 0.8 and Qm = 0.3. The solid line 

mark the upper limits when we assume T = 8.5M0g54 (Nevalaincn 

et al. 2000). If the assumed T - M relation is T = S.OM^f^ 
(Evrard ct al. 1996) the upper limits on Lo change (dotted line). 
We also show the case when a pre-heating like model with two 
different exponents is assumed (dashed line), Lx = Lo * T°' for 
T > 1 keV and Lx=Lo* Tl^ for T < 1 keV. We assumed /3 = 2a 
and the Nevalainen et al. (2000) T — M relation. This shows how 
the constraints are not very sensitive to clusters with T < 1 keV. 



as. If we parametrise the luminosity as La; = LoT", then 
we can constrain Lo as a function of a by just requiring 

that the cluster power spectrum stays below the measured 
power spectrum (see figure 16). Although we have not dis- 
cussed much about the T — M relation, it will also play a 
role in the definition of the upper limit. To get the temper- 
ature in the Lx — T relation we need to assume some form 
for T — M relation. Recent observations (Nevalainen et al. 
2000), seem to indicate that this relation differs slightly from 
the self similar scaling relation. However, as in the case of 
the Lx — T, there is some scattering on the T — M relation. 
We have included this uncertainty into our calculations by 
using two different T — M relations. In the first case we use 
the fit to the observed T and M obtained by Nevalainen et 
al. (2000), T = 8.5Mi°5^* keV and in the second case we use 
a self-similar form (Evrard et al. 1996), T = S.OMj/^ keV. 
The results on Lo and a are shown in figure 16. 

When we compare the models on the two upper limit 
curves with the observed L^—T relation (figure 17), we find 
that the models on the upper limits are marginally consis- 
tent with the observations. The model with the self-similar 
T — M relation produces a constrain which agrees better 
with the observations. 



Figure 17. L^—T relations for the models lying in the two upper 
limits of figure 16. The dotted line corresponds to the models in 
the solid line of figure 16 while the dashed lines correspond to the 
models in the dotted line curve in figure 16. The solid line shows 
the prehetaing model of Babul et al. (2002). 



7 DISCUSSION 

We have found that in order to explain most of the con- 
straints obtained in as and f2m, we need a rather steep 
Lx — T relation. We have changed these relations within 
the observed limits (and even more, in the case of the pa- 
rameter a) and we found that the as = 0.8 Q,m = 0.3 model 
can be accommodated only in the most extreme cases. Using 
the pre-heating model of Babul et al. (2002), our upper limit 
excludes most of the current constraints in ag for Sl,„ ~ 0.3. 
If as > 0.8 for Qrn = 0.3, we may say that there is an incon- 
sistency between the above preheating model and the upper 
limit imposed by the RASS diffuse background. 

As noted by Voit et al. (2002), the observed Lx—T rela- 
tion changes when the luminosities are corrected by cooling 
fiow effects. Our results suggest that cooling flows can be 
important (specially at low redshift which is the range dom- 
inating the power spectrum in this scale). A L^ — T which 
is obtained without correcting for cooling flows will be af- 
fected by a larger scatter due to the presence of the clusters 
with cooling flows (they will increase the average luminos- 
ity of clusters and therefore the value of Lo). If cooling flow 
clusters are not representative of the average population of 
clusters, they should be excluded from any fit to the Lx—T 
relation or corrected their luminosities (and average tem- 
peratures) by the effect of the cooling flow in order to get a 
representative — T relation (see Fabian et al. 1994 for a 
discussion on the effect of cooling flows in the L^ — T rela^ 
tion). 

In figure 18 we present one of the extreme cases for 
which the model as = 0.8, Qm = 0.3 is just below the up- 
per limit (solid line). This model corresponds to the values 
Lo = 2.3 X 10''^/i~^ ergs~^ and a = 3 and is in the dotted 
curve of figure 16 (i.e. when we assume the self-similar re- 
lation T — S.OM^g^ keV). If we compare this extreme case 
with real observations of the Lx — T relation, we see that this 
model still produces a good fit to the data (but we need to as- 
sume that T oc Mjg^). All the models lying on the right hand 
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side of the solid line overpredict the observed RASS power 
spectrum for as ~ 0.8, Qm ~ 0.3. The pre-heating model of 
Babul et al. (2002) (bottom dotted line) is marginally consis- 
tent with the upper limit above 4 keV but exceeds the upper 
limit below that temperature. On the other hand, the pro- 
heated model of Voit & Bryan (2001) is still compatible with 
the upper limit but is well within the limits. There is little 
room for additional changes in the Lx — T and T — M relar 
tion without coming into serious contradictions with the ob- 
served relations. As a consequence, our results suggest that 
it is very difficult to accommodate high as models {as > 0.8) 
within our upper limit. Recently, a combined analysis of 2dF 
plus CMB data has produced now constraints on as which 
suggest that as ~ 0.7 for Qmh = 0.21 (no tensor model) and 
as ~ 0.6 for ilmh — 0.16 (tensor model) (see Efstathiou et 
al. 2002). Other recent works also suggest a low value for 
0-8 (e.g Allen et al. 2002, Seljak 2002, Bahcall et al. 2002). 
These values of as are in very good agroomont with our up- 
per limits. If the value of as is significantly larger than 0.8, 
(for = 0.3) then, we can conclude that there is a serious 
contradiction between the X-ray data and the model. If we 
adopt the realistic position that the inconsistency must be 
in the model, then, wo should revise our assumptions in or- 
der to resolve this puzzling situation. 

The only assumption we have not checked in our model is 
that about the mass function. As mentioned in the introduc- 
tion, the modelling of the cluster power spectrum involves 
the whole population of clusters. This fact makes the cluster 
power spectrum a very sensitive probe of the cluster mass 
function. 

This opens the exciting possibility of (assuming the cosmo- 
logical model is known) constraining the cluster mass func- 
tion once wc know accurately the Lx — T and T — M relations 
and a good estimation of the X-ray background power spec- 
trum is available. 

There is an important point which has not been ex- 
plored in this paper but that wc plan to study in a sub- 
sequent work. It is the scattering of the — T relation. 
Real clusters will show an intrinsic scattering in this rela- 
tion (clusters with the same temperature may have different 
luminosities around the average Lx)- Since the power spec- 
trum of clusters is proportional to the square of the multi- 
pole decomposition (equations. 1 and 3), those clusters with 
the same T but Lx larger than the mean will contribute 
more to the power spectrum than the clusters with Lx lower 
than the mean (but the same T). That is, the average effect 
of the population of clusters with the same T (but differ- 
ent Lx) is not the same as the effect of the population of 
clusters with the same T and average luminosity Lx (due 
to the quadratic dependence of Ci). This fact may change 
the results presented in this work and need to be studied in 
more detail. 



8 CONCLUSION 

We have shown that the use of the measured X-ray power 
spectrum as an upper limit produces strong constraints on 
the cosmological model. The constraints have the advan- 
tage of being very robust since they were obtained under 
a minimum number of assumptions. The uncertainty in the 
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Figure 18. Scaling model for the curve in the upper limit 
Lo = 2.310*^ ?i~^erg/s and a = 3 (solid line). The model is com- 
pared with a recent compilation of temperatures and luminosi- 
ties (Novicki et al. 2002) at low (dots) and intermediate redshift 

(open triangles) . The two dotted lines show the preheating models 
of Voit & Bryan (2001) (top) and Babul et al. (2002) (bottom). 
Note the large scattering in the data. 



Figure 19. Upper limit on erg — f2m compared with a recent con- 
straint in these parameters from CMB (Melchiorri et al. 2002). 
The X-ray backgound upper limit excludes half of the CMB con- 
straints. 



assumptions cam be easily incorporated into the model to 
produce optimistic and pessimistic constraints. Approaches 
like the one presented in this work arc needed to comple- 
ment the constraints obtained using other methods (see fig- 
ure 19). To obtain the upper limits on as — Q,m presented in 
this work we do not even need to detect the clusters. The 
global contribution of all the clusters (above or below the So- 
level) contributes to the power spectrum although we have 
seen that the cluster population with intermediate masses 
and z < 0.2 dominates the power spectrum in the range of 
scales considered in this work. 

In this work, we have not considered the clustering contri- 
bution (two-halo). This contribution may be of the order of 
20% — 30% at the scales relevant for this work (Lahav et 
al. 1997, Komatsu & Kitayama 1999). The upper limits pre- 
sented hero will be even more significant if this contribution 
is incorporated into the analyses. 

Due to the point source contamination and large pixel size, 
the constraints on as — 0.m arise in the range / < 100 — 200. 
However, this fact has an advantage. At low Vs, the cluster 
power spectrum is independent of the electron density pro- 
file (or the internal structure of the cluster) and only the 
total luminosity of the clusters matters for modelling the 
power spectrum. This reduces dramatically the number of 
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assumptions made in the model. One may use this fact to 
make studies of the physics of the intracluster plasma (like 
pre-heating) which could constrain the Lx —T relation. 
We found that, in order to accommodate cosmological mod- 
els like (Tg = 0.8 Q,m = 0.3 wo need to tune the T — M and 
the Lx — M relation to the limits of the observational con- 
straints. This fact suggest that lower as models are favoured 
by the upper limit constraint set by the X-ray background. 

The approach presented in this work can be extremely 
useful with higher resolution and better signal to noise data 
{Chandra, and particularly XMM-Newton) . Although the 
data in the ROSAT R6 band has the brightest point sources 
removed, the power spectrum still contains some contribu- 
tion coming from the non-removed point sources and their 
clustering. Future data from Chandra and specially XMM- 
Newton will allow a better point source subtraction. The 
constraints obtained from these new data will be much more 
significant since one can probe smaller scales with a sig- 
nificantly smaller contamination from point sources. In the 
small scale range, the cluster power spectrum becomes even 
more sensitive to the physics of the plasma (angular exten- 
sion of the emission and evolution of the L^ — T relation) . In 
particular, we have seen that the shape of the power spec- 
trum at small scales (large multipoles) is better described by 
density profiles having high dumpiness factors (p « 10—20). 
A low-contaminated point source power spectrum could lead 
to dramatic constraints on the Lx — T relation. These con- 
straints can have serious implications for the physics of the 
intracluster plasma (pre-heating, cooling fiows). These facts 
make worth extending the analyses presented here to smaller 
scales with a high sensitivity where the point source removal 
can be performed to very low fluxes. 
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